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ABSTRACT 

Constraints on cosmology from recent cosmic shear observations are becoming increas- 
ingly sophisticated in their treatment of potential systematic effects. Here we present 
cosmological constraints which include modelling of intrinsic alignments. We demon- 
strate how the results are changed for three different intrinsic alignment models, and 
for two different models of the cosmic shear galaxy population. We find that intrinsic 
alignments can either reduce or increase measurements of the fluctuation amplitude 
parameter erg depending on these decisions, and depending on the cosmic shear survey 
properties. This is due to the interplay between the two types of intrinsic alignment, 
II and GI. It has been shown that future surveys must make a careful treatment of 
intrinsic alignments to avoid significant biases, and that simultaneous constraints from 
shear-shear and shear-position correlation functions can mitigate the effects. For the 
first time we here combine constraints from cosmic shear surveys (shear-shear correla- 
tions) with those from "GI" intrinsic alignment data sets (shear-position correlations). 
We produce updated constraints on cosmology marginalised over two free parameters 
in the halo model for intrinsic alignments. We find that the additional freedom is well 
compensated by the additional information, in that the constraints are very similar 
indeed to those obtained when intrinsic alignments are ignored, both in terms of best 
fit values and uncertainties. 

Key words: cosmology: observations - gravitational lensing ~ large-scale structure 
of Universe - galaxies: evolution 



1 INTRODUCTION 

Weak gravitational lensing is one of the most promising ob- 
servational tools available to cosmologists for studying the 
recent accelerated expansion of the Universe. Images of dis- 
tant galaxies appear distorted due to the bending of light 
as it passes through the gravitational potential of the inter- 
vening matter. The image distortion can be described as a 
shearing of the original galaxy shape. This "cosmic shear" 
can be exploited to study the matter distribution in the 
Universe and the growth of structure. Cosmic shear has the 
potential to be one of the mo st powerful probes both of 
dark matter and of dark energy (|Peacock fc Schneideij[200^ : 
lAlbrecht et al] bood ') which together are thought to make 
up about 95% of the energy budget of the Universe (e.g. 
ISpergel et al]|2007l '). 

For any individual galaxy it is impossible to separate 
the small cosmic shear distortion from the intrinsic elliptic- 
ity of the galaxy shape. However, the light from physically 
close galaxies will follow a similar trajectory to the observer. 



passing through spacetime curved by the same gravitational 
fields. Therefore these galaxies will acquire the same cos- 
mic shear distortion. If it is assumed that galaxy shape and 
orientation are randomly assigned across the sky then the 
cosmic shear signal can be retrieved by averaging the ellip- 
ticities of a number of galaxies close on the sky. 

In practice this assumption of randomly distributed 
galaxy shapes is unrealistic. Spatially localised galaxies are 
expected to have formed within the same large-scale gravita- 
tional field, which i s likely to cause an alignment in their in- 
trins i c ellipticities (iHeavens fc PeacockllipS^: Catelan et al.l 
I2OOII : |jinell2002l : lAubert et al.ll2004 ''). See'S chafeij (|2008h for 
a recent review. This intrinsic alignment (lA) effect will 
appear as a systematic error in estimates of cosmological 
parameters extracted from cosmic shear data, unless accu- 
rately taken into account. 

Two types of lA affect the measured cosmic shear 
signal. Physically close galaxies form in the same large- 
scale gravitational potential so share a preferred ellipticity 
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orientation. This is the Intrinsic-Intrinsic (II) correlation. 
Background galaxies will be lensed by foreground gravita- 
tional potentials which govern the orientation of foreground 
galaxies. This causes an anti-correlation between the fore- 
ground/background galaxies known as the Gravitational- 
Intrinsic (GI) correlation. 

Broadly, there are two approaches to dealing with lAs. 
The first, known as "nulling", places steps in the weak lens- 
ing pipeline which remove the lA signal. The effect of II cor- 
relations can be removed by downweighting physically close 
galaxy pairs jHevmans et al.l 20oj: King fc Schneiderll2003l : 



iHevmans fc Heavensll2003l : iTakada fc WhitdbOoj TThe GI 
term is more problematic as it affects all pairs of galaxies 
which are not physically close. In principle a particular lin- 
ear combination of tomographic shear power spectra can be 
used to remov e the GI signal if the redshifts o f the galax- 
ies are known (| Joachimi fc SchneideilboOSl . |2009| ) . The other 
approach tries to model the lA signal we would expect for 
a particular survey. These lA contributions can then be in- 
corporated into the predictions for the measured shear sig- 
nal and any free param eters in the lA model can be varied 
and m arginalised over (Alb recht et~ar 2006; B ridle fc Kin3 
I2OO7I : [Bernstein 2009 : Joachimi fc Brid le 2009,). This is the 
approach we investigate. To carry this out it is neces- 
sary to make physically motivated models for I As. Ignor- 
ing lAs c ompletely will bias est i mates of cosmological pa- 
rameters (iHir ata fc Scliak" 2004; M andelbaum et al.l l2006l : 
iHirata et allBoO?. : Bridle fc King. 200"7t ). In addition, other 
cosmological data sets, such as galaxy surveys, can be used 
to gain empirical knowledge about the I A effect, as well as 
its variation with luminosity, colour and any other variables. 

Some of the modelling approaches make simultaneous 
use of the galaxy shear - galaxy shear correlation function 
and the galaxy position - galaxy shear corre lation function, 
as measured from the same imaging survey (* Albrec ht et al.l 
[2OO& : Berns tein 2009; Joachimi fc Bridle 2009). Where 
"shear" has contributions from both gravitational and in- 
trinsic shear. This was first identified as an important ad- 
ditional st a tistic for learning about galaxy formation by 
IHu fc JainI l|2004 ) and was rec ently sugg ested as a power- 
ful tool for remo ving lAs by .Zhang et al.l (12009 '). which was 
demonstrated in I Joachimi fc Bridld ( 2009h . 

Cosmic shear was first detected observationa,lly a 
decade ago jBacon et al. _ 2000l: Van Waerbeke et al] I2OO0I : 



IWittman et al.ll2000l : iKaiser et al.ll200rf ) and the strength of 
the signal as a function of galaxy separatio n has since been 
mea sured by many team s, most re cently by Fu et al.l (12008 



mea sured by many team s, most re centlv bvli^u et al 
and iMassev et al.l ((20o3) (see also ISchrabback et al 



A compilation of recent cosmic shear results including a ho- 
mogeneou s additional treatment of systematics was made 
public by iBeniamin et al.l ()2007l ) . The use of photometric 
redshift information is now starting to allow the signal to 
be ev aluated as a funct ion of galaxy redshift (jMassev et al.l 
12007; Schra bback et alJ , 2009). The cosmological constraints 
in these recent papers are calculated assuming there are no 
I As. 

However, if lAs are non-negligible then ignoring them 
will lead to a bias on cosmological parameters. The bias 
in (78 was estimated to be be tween 1 and 20 p e r cen t 
for a CFHTLS-hke survev bv fMandelbaum et all (l2006l l 



approximate Fisher matrix prediction. iBridle fc KiiiS (|2007l ) 
showed that for future surveys the equation of state param- 
eter w may be biased by order unity if lAs are significant 
and ignored. 

iHevmans et al.1 (|2004 ) used COMBO-17 measurements 
of the II signal from shear-shear correlations between galax- 
ies close in redshift jointly with cosmic shear RCS and 
VIRMOS-DESCART shear-shear correlation data to con- 
strain a simple intrinsic alignment model amplitude simul- 
taneously with cosmological parameters. They marginalised 
over possible amplitudes of the II signal and found the fitted 
value of (78 was reduced by 0.03 relativ e to the v a lue fo und 
when intrinsic alignments were ignored. IFu et all (|2008l ) es- 
timated the amplitude of the GI signal marginalised over 
a range of as, Qm values allowed by the cosmic shear data. 
They found it to be consistent with zero, and this conclusion 
held when various different scale ranges were used for each 
of the fits for cosmology or for the lA amplitude parameter. 
They con cluded that they found n o evidence for a non-zero 
GI signal. ISchrabback et al] (|2009l ) took significant steps to 
reduce the contamination of their cosmic shear constraints 
by lAs by excluding the auto-correlations for the 5 narrow- 
est redshift bins and excluding luminous red galaxies (LRGs) 
from their analysis. They also showed results using only the 
autocorrelations, which were similar to those excluding the 
autocorrelations results, leading them to conclude that in- 
trinsic alignments are not a significant contaminant. 

Many current and future surveys are planned with 
cosmic shear as a major design driver, in particular from 
the ground the Canada-France-Hawaii Telescope Legacy 
burvey (CFHTLSfl the Kilo-Degree Survey (KIDS), the 
Panoramic Survey Telescope and Rapid Response Sys- 
tem (Pan-STARRSfl surveys, the Dark Energy Survey 
(DESfl the Large Synoptic Survey Telescope (LSSTf], and 
space missions Euclicfl and the Joint Dark Energy Mission 
(JDEMfl. 

The majority of observational constrain ts on galaxy 
lAs h ave been carried out at low redshift (iBrown et al.l 
l2002l : IHevmans et al. 1 12004 iMandelbaum et al.l I2OO6I ) but 
these have recently been extended towards the red- 
shifts relevant to cosmic shear surveys (jHirata et akllioOTl : 
IMandelbaum et al] |2009| ). These studies have used spec- 
troscopic galaxy surveys to calculate the galaxy position - 
galaxy shear correlation function and/or the galaxy shear 
- galaxy shear correlation function. Constraints calculated 
from these correlation functions have been fitted with sim- 
ple models for lAs using a fixed cosmological model. 

In this paper we calculate constraints on cosmological 
parameters from cosmic shear data, taking into account the 
likely contamination from lAs for a range of lA models. We 
simultaneously fit an lA model and cosmological parameters 
to lA correlation functions. We then perform a joint analysis 
of cosmic shear and lA data, varying both cosmology and 
parameters within the lA model. 

The paper is organised as follows; Section [2] reviews the 



e y - . 

and IHirata et all (l2007h ISchneider fc Bridle! (|201Ci ) and 
IMandelbaum et al] ( 20091 ) obtained a similar result in an 
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cosmic shear formalism and some basic models of lAs; Sec- 
tion |3] presents constraints from the shear-shear correlation 
data; Section |3] introduces an lA parameterisation and con- 
straints from shear-position data; Section [5] uses shear-shear 
and shear-position data to produce joint constraints on lA 
and cosmological parameters and we conclude in Section [6] 
Throughout we assume a fiat, ACDM cosmology. When 
we refer to fiducial values for cosmological parameters we 
set the present day amplitude of linear mass fluctuations 
(T8= 0.751, the baryon density Qh~ 0.05, the matter density 
i^m= 0.3, and the Hubble parameter defined by Ho ~ lOO/i 
km s~^ Mpc~^ as h= 0.7. We assume a Ifarrison-Zel'dovich 
primordial power spectrum slope Us — 1. 



2 COSMIC SHEAR AND INTRINSIC 
ALIGNMENTS 

In this section we review the cosmic shear formalism before 
describing the origin of lAs and their effect on the mea- 
sured cosmic shear signal. We then summarise three ba- 
sic approaches to modelling lAs, beginning with the linear 
alignment model, then its non-linear extension and finally 
the halo model of lSchneider fc Bridle! (|2010l '). 



2.1 Cosmic Shear 

Cosmic shear probes the expansion history and the growth 
of structure in the Universe using weak gravitational lensing. 
The effect is rather small (~ 1%) for most galaxies especially 
when compared to the intrinsic ellipticity of a single galaxy 
(~ 20%). Therefore it is necessary to average over many 
galaxies to detect a coherent signal. A convenient statistic 
is the Fourier transform of the two-point shear correlation 
function (^r) between pairs of galaxies, the cosmic shear 
power spectrum Cp''. It is well approximated by projecting 
the three-dimensional matter power spectrum Ps{k;x) ^ 
redshift corresponding to a line-of-sight comoving distance 
X, weighted by the lensing kernel, W{x) 



(1) 



where in a flat universe 



W{x) = l^^Hla~\x)x r" dx. n{xs)^^^-^ (2) 



and a is the scale factor, xh is the comoving distance to 
the Hubble horizon, and the function n{xs) is the selection 
function of source galaxies per unit comoving distance, nor- 
malized as / dxn{x) = 1- 

The two-point shear correlation function is determined 
by the cosmic shear power spectrum via 



(3) 



^E{e) = — / AiiCiMW), 



where Jo is the zeroth order Bessel function of the first kind. 



2.2 Intrinsic Alignments 

If galaxies were randomly oriented on the sky then the cos- 
mic shear signal would be relatively easy to observe because 



the random intrinsic contributions to the observed elliptici- 
ties would average away to zero. However this is not expected 
to be the case. Galaxies form within large-scale gravitational 
potentials, producing some "intrinsic alignment" between 
galaxy shear. There are two ways in which lAs contribute 
to the observed shear power spectrum: 

The II correlation: Intrinsic-Intrinsic galaxy alignments 
occur because physically close galaxies form in the same 
tidal field, causing an alignment of their halos or angular 
momentum vectors. As a result the galaxies point in the 
same direction. This alignment produces an increase in the 
measured shear power spectrum. 

The GI correlation: Gravitational-Intrinsic alignments 
are a cross term between intrinsic ellipticity and cosmic 
shear. The intrinsic ellipticity of a galaxy is aligned with 
the density field in which it forms and this field in turn con- 
tributes to the lensing distortion of more distant galaxies. 
This double role causes an anti-correlation between galaxy 
ellipticities, because the closer galaxy may be expected to 
point towards the local overdensity while the distant galaxy 
is stretched tangentially around the overdensity. This leads 
to a suppression of the measured power spectrum. 

The measured shear power spectrum C; therefore arises 
from genuine cosmic shear, combined with I A terms pro- 
duced by the II and GI correlations 



(4) 



The strength of the II term depends strongly on the 
depth of the survey. A deep survey may contain galaxies at 
many different distances from the observer, so two galaxies 
which appear close on the sky are most likely to be far apart 
along the line of sight and therefore not physically close. 
Since the opposite is true for a shallower survey then the lA 
is stronger. In addition the lensing strength is smaller when 
all distances are reduced, so the cosmic shear contribution 
is small for a shallow survey. 

The II and GI lensing power spectra are calculated anal- 
ogously to Cp^ 



Jo 

jEE 



n\x) 



P^i''{k,x)Ax 



2W{x)n{,x) 



(5) 



(6) 



where k = l/x, P^i {k,x) P^^/(fc,x) are the intrinsic- 
intrinsic and gravitational-intrinsic power spectra respec- 
tively. 

In this paper we illustrate the effect of lAs on cosmo- 
logical constraints usin g data from the 100 Square Degree 
Weak Lensing Survey ( Benjamin et al.ll2007l ). We compare 
constraints from four different approaches to lAs: 

(i) Ignoring lA effects 

(ii) The linear alignment model of lAs 

(iii) The non-linear alignment (NLA) model of lAs 

(iv) The Halo Model of I As. 

lAs can be studied with probes other than the usual 
cosmic shear measurement of shear-shear correlation func- 
tions. Section[4]below we use data from the Sloan Digital Sky 
Survey (SDSS) galaxy redshift surv ey shear-position corre- 
lation functions (|Hirata et al.|[2007l ) to constrain the ampli- 
tude, scale and luminosity dependence of lAs. These results 
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can be combined with models of lA in cosmic shear pre- 
dictions to produce less biased constraints on cosmological 
parameters. Details of the lA parameterisations are given 
below. 

2.2.1 Linear Alignment Model 

A popular simple method for modelling lAs is known as the 
Linear Alignment (LA) model because it assumes that the 
intrinsic ellipticity of galaxies is proportional to the curva- 
ture of the primordia.1 larg e scale potential jCatelan et al.l 
l200ll : lHirata fc Seliakl l2004li. This is thought to be most rel- 
evant for elliptical galaxies. 

In the linear alignment model the II and GI power spec- 
tra are found to first order to be 

= ^Pt{k) (7) 

Ps,,i{k) = ~^P^'"{k) (8) 

where Pg^"{k) is the linear theory matter power spectrum, 
Ci is a normalization constant, D{z) = (1 -I- z)D{z) where 
D{z) is the growth factor normalized to unity at the present 
day, and p{z) is the m ean matter density of the Universe as 
a function of redshift l|Hirata fc Seliaklliooj '). 



radius vector of the parent halo. Shear two-point correlation 
functions then consist of two components: a "one-halo term" 
arising from pairs of galaxies within the same halo, and a 
"two-halo term" due to pairs in different halos. 

For lAs the one-halo term comes from the radial align- 
ment of the satellite galaxies and the two-halo term is domi- 
nated by the correlation of the central galaxies in each halo. 
Within this model one- and two-halo terms appear in both 
the II and GI lA power spectra 

pEE _ pEE,lh . pEE,2h ,g\ 

Ps,^: ^Pt^i +P!I;i. (10) 

In the case of the II power spectrum, the two-halo term 
can be further broken up into terms due to the correlation of 
satellites with satellites, satellites with centrals and centrals 
with centrals, and, in the GI case, the two-halo power spec- 
trum can be separated into a central and a satellite term. 
However, it was found (jSchneider fc Bridle 2010 ) that the 
two halo terms are dominated by the central with central 
(II) and central (GI) terms, allowing us to ignore two-halo 
terms involving satellites. This greatly simplifies the equa- 
tions, allowing us to write down the two-halo terms 

pEE,2H ^ ^pU^^k) (11) 



2.2.2 Non-Linear Alignment Model 

This linear alignment model is likely to hold on large scales 
but it takes no account of non-linear growth of structure. 
One somewhat ad-h oc method that has been employed 
l|Bridle fc Kind [20071 ) to rectify this deficiency is to substi- 
tute the non-linear matter power spectrum P^'(fc), instead 
of the linear theory matter power spe ctrum Pj'"(fc ) , into 
Eq.[71and Eq.[Hl In this work we use the lSmith et al.l (|2003l ') 
prescription for the non-linear matter power spectrum. We 
call this approach the non-linear alignment (NLA) model. 

2.2.3 Halo Model 

The halo model of galaxy clus tering 

IScherrer fc Bertschingeij Il99ll : IScoccimarro et all 1200 ll ) 
is a simple but effective predictor for galaxy clustering 
statistics. To first order in this model the Universe consists 
of dark matter halos, distributed according to linear theory. 
Each halo is described by a mass, drawn from a mass 
function, and a density profile. The forms of these functions 
are gen erally taken fro m ave rages of n-body simulations 
(see e.g. ICoorav fc SIiethH20()3 . for a review). 

Applying the halo model to lAs provides a physically 
motivated way to model small-scale fea tures of the lA sig- 
nal. A simple example was provided in ISchneider fc Bridle! 
(|2010l ) which we describe and use here. Each dark matter 
halo is taken to have a single central galaxy, positioned ex- 
actly at the centre of the spherical halo mass distribution. 
These central galaxies are assumed to have ellipticities de- 
termined by the large-scale density perturbations according 
to the linear alignment model. The satellite galaxy number 
density follows the density profile of the halo and their ellip- 
ticities are defined according to a distribution in the angle 
between the satellite major axis and the three-dimensional 



p,'^. = -^p'r{k) (12) 

i.e. the II and GI two-halo terms are just the same as the 
full linear alignment model of lAs. 

[Schncidcr fc Bridle (2010,) provide fitting formulae for 
the one-halo power spectra 

pEE,lh,,. _ 2 (k/pi)^ , . 

P^,¥(k) = ~'ysc.^.YT{kj^ (14) 

where yscaie is the degree of alignment and determines 
the amplitude of the Ih power spectra, that could be 
set by comparison with simulations, and for the functions 
Pi( z) (for i = 1,2,3) we use the fitting formulae provided 
in ISchneider fc Bridle! ()20ld ). that are based on a halo 
model wi th redshift dependence determined by the Sheth- 
Tormen (jSheth fc Tormenlll999l ) model for the halo mass 
function and bias as well as the linear growth function in 
their fiducial cosmology. 

The lA power spectra and Ps,ji enter both the 

shear-shear and shear-position calculations. In principle 
they depend on cosmology through the contents of the one- 
and two-halo terms, given in equations 1111 to [T4l The two- 
halo term contains the linear theory matter power spectrum 
and the growth rate, both of which depend on the values of 
cosmological parameters. The one-halo term is derived from 
a consideration of the number halos as a function of halo 
mass, and a match to the redshift distribut ion o f the s urvey 
in question. It was shown in Schn eider fc Bridlel (|2010| ) that 
the dependence on the survey redshift distribution is weak 
and thus we neglect it here. The dependence of both terms 
on cosmology could be included, and this would make sense 
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if we had complete faith in the models. We have the alter- 
native of considering the lA power spectra for our fiducial 
model as a template representing a reasonable but poorly 
understood guess at the contribution fro m lAs. For this work 
we us e the one-halo term as given in [Schneider fc Bridle! 
l|2010l ) and ignore any cosmology dependence within it. We 
choose to fix (Tg within the two-halo calculation to the fidu- 
cial value. Therefore the relative amplitude of the one- and 
two-halo terms stay roughly constant and the overall am- 
plitude of the sum is varied (see below). In this paper we 
choose to vary the matter density inside the linear theory 
matter power spectra used in the two-halo terms in equa- 
tions[TT]and[T2]and in the Sim contribution to the normalised 
growth in the denominator of these equations. 

To include the colour and luminosity dependence we 
generalise our halo model equations [9l and 1101 



f r,EE,lh . pEE,2h 



fr 



fr 



(15) 



(16) 



which scales the amplitude of both one- and two-halo terms 
by the same factor, given in square brackets. The factor in 
square brackets is squared in the first equation and not in 
the second equation, to mimic a simultaneous modulation of 
the one-halo scaling parameter 7scaio and the two-halo am- 
plitude parameter Ci . Throu ghout we retain the fiducia l val- 
ues of 7acaio = 0.21 following ISchneider fc Bridl3 (|2010|') and 
Ci = 5 X lG-^*(h^Mf^Mv,c^)-^ following iBridle fc Kind 
l|2007t ). 

L/ Lq is the normalised luminosity of the data bin and 
fr is the fraction of red galaxies in the data bin. This power 
law i n luminosity is equiv alent to that used in the power law 
fits of lHirata et all \2m% . The motivation for multiplication 
by the red fraction fr is that we assume only red galaxies 
have lAs. These equations could be generalised to have dif- 
ferent lA amplitudes for red and blue galaxies by adding a 
term proportional to the blue fraction {1 — fr) with a differ- 
ent variable amplitude parameter. However, we defer such 
modelling to future work. A and /3 are free parameters. Note 
that the above equations reduce to the basic halo model for 

A = fr = 1,^ = 0. 



3 SHEAR-SHEAR CORRELATIONS 

In this Section we summarise the cosmic shear data we use 
and compare it with predicted correlation functions with the 
various lA models. We then make a first calculation of the 
impact of the models on constraints on the amplitude of 
matter clustering and the matter density of the Universe. 
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Figure 1. The upper panel shows galaxy redshift distribution, 
n(z), for the four surveys comprising the 100 deg^ cosmic shear 
dataset. The lower panel shows the n(z) difference with respect to 
CFHLTS, n(z)-n(z)cFHTLS- CFHTLS is the solid line, GaBoDS 
is dot-dashed, RCS is dotted and VIRMOS-DESCART is dashed. 



covers an area of 22 deg^, reaching a depth of i' = 24.5 
(jHoekstra et al.ll2006l) . There is 13 deg^ of data from the 
GaBoDS survey which uses objects which lie in the interval 
R G [21.5,24.5] (|Hetterscheidt et al.|[2007l ). The RCS data 
covers 53 deg^ with a limit ing magnitude of 25.2 in the Rc 
band (|Hoekstra et al.ll2002t ). The VIRMOS-DESCART data 
has an effective a rea of 8.5 deg^ an d a limiting magnitude 
of Iab = 24.5 (Le Fevre et al.ll2004h . 

Throughout this paper we use t he redshift distributions 
for each data set as given by Eq. 9 of [Benjamin et al.l (|2007l ) 
using the parameters given in the upper section of their Ta- 
ble 2 for the high confidence regime (fit ted to photometri c 
redshifts in the range 0.2 < z < 1.5 from lllbert et ahlbood ). 
iBeni amin et al ] (|2007h took into account uncertainties in the 
photometric redshift distributions by sampling from tuples 
of these parameter values and selecting the best fit values for 
each of the three parameters. For computational practical- 
ity we here ignore uncertainties in the photometric redshift 
distributions. We therefore next describe in more detail the 
relative differences between the redshift distributions of the 
different surveys, ignoring the uncertainties. 

Since the impact of lAs on each dataset hinges so much 
on the redshift distribution we show the redshift distri- 
butions and the difference in distribution, with respect to 
CFHTLS, for each dataset in Fig.[Tl We see that VIRMOS- 
DESCART and GaBoDS have the greatest number of low 
redshift galaxies and VIRMOS-DESCART shows a signifi- 
gant excess at z ~ 1.4. 



3.1 Cosmic Shear Data 

The dataset used to constrain models in this paper is the 10 
Square Degree weak lensing survey (jBeniamin et al.ll2007h . 
This combines data from the Canada-France-Hawaii Tele- 
scope Legacy Survey (CFHTLS)-Wide, the Garching-Bonn 
Deep Survey (GaBoDS), the Red-sequence Cluster Survey 
(RCS) and the VIRMOS-DESCART surveys. 

The CFHTLS- Wide data included in this compilation 



3.2 Correlation Functions Including Intrinsic 
Alignments 

Fig.[5]shows the correlation functions of the four constituent 
surveys of the 100 deg^ dataset. We overlay the predicted 
correlation functions for each of the lA models we consider. 
The correlation functions are calculated for the fiducial cos- 
mology, with no at tempt to find th e best fit to each survey. 
Note that Fig. 1 in lBeniamin et al] (|2007l ) plots similar cal- 
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Figure 2. Shear correlation function, g^i data from the four surveys comprising the 100 deg^ cosmic shear dataset along with correlation 
functions for each of the four lA models using the fiducial cosmology and the appropriate n(z) for each survey. The correlation function 
without lAs is shown by the solid line, the linear alignment model as dotted line, the NLA model as dot-dashed and the halo model as 
the dashed line. The data, with 68% errors, appears as crosses. 



culated correlation functions but uses a different, best-fit, 
as for each survey. 

Overall the GI contribution dominates the lA signal 
and has the effect of reducing the shear-shear correlation 
signal; hence all the correlation functions which include an 
lA component are mostly lower than the equivalent corre- 
lation function without lAs (solid line). As expected, the 
halo model (dashed line) agrees with the LA (dotted line) 
and NLA (dot-dashed line) models at large angular separa- 
tion where the two-halo term and linear theory dominate. 
The halo and LA models diverge at small angular separation 
when the additional one-halo term becomes important and 
the shear-shear correlation amplitude is reduced due to the 
negative contributions from the GI term. The NLA model is 
intermediate between the LA and halo model at small angu- 
lar separations. The increase in the halo model correlation 
function at intermediate scales is caused by an increased II 
contribution on these scales. In the halo model of lAs the 
galaxies align towards the center of dark matter halos, and 
as a result are also aligned with each other. Therefore the 
one-halo term for the GI signal dominates on a scale corre- 
sponding to the typical distance of galaxies from the halo 
center, whereas the one-halo term for the II signal occurs at 
the typical halo diameter. This explains why the GI term 
dominates at small scales and the II term has a bigger im- 
pact at intermediate scales. 

We now compare the predicted correlation functions 
for the different surveys. The only difference between the 
predictions comes from the different redshift distributions 
used. The lA signals dominate over the cosmic shear signal 
at low redshift. A large number of physically close galax- 
ies produces a large II effect. Surveys with a large num- 
ber of low redshift galaxies will have a large proportion 
of these physically close strongly aligned galaxies. Fig. [T] 



shows that VIRMOS-DESCART (dashed line) and GaBoDs 
(dot-dashed line) have the largest proportion of low redshift 
galaxies, so we expect the strong increase at intermediate 
scales in Fig. [S] Conversely CFHTLS has the lowest pro- 
portion of low-z galaxies and we see little increase in its 
correlation function due to an II contribution. 

3.3 Impact on erg constraints 

We are interested in the biasing of cosmological constraints 
when lAs are erroneously ignored. We compare predicted 
and observed correlation functions and produce values at 
a range of amplitudes of matter clustering, parameterised 
by (jg. These are used to calculate a best fit erg and cor- 
responding 68% confidence limits. This was done for each 
of the four cosmic shear datasets separately and the anal- 
ysis was repeated using the four different I A approaches. 
The results are shown in Fig. |3] The right-hand panel 
shows the 68% confidence limits using the combined 100 
deg^ data, obtained by multiplying together the probabili- 
ties Pr((jg) (X exp(— x'^/S) from each survey. In this analysis, 
all other cosmological parameters are held at their fiducial 
values and we fix A = 1, /3 = and /_b = 1. In all analy- 
sis quoted in this paper, we kept erg constant, at its fiducial 
value, within the lA models themselves which can be scaled 
by an overall amplitude parameter A. As usual, the cosmic 
shear correlation function has a strong dependence on ag. 

The results show a dependence of the measured erg on 
the lA model. Ignoring lAs (horizontal/ vertical crosses) pro- 
duces results which are consistently below the results when 
lAs are taken into account, irrespective of the I A model used. 
This behaviour is consistent with the correlation functions 
seen in Fig. (2] The dominance of the GI term has the effect 
of lowering the overall shear-shear correlation function pre- 
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Figure 3. 68% confidence limits on the amplitude of fluctua- 
tions as for each survey within the 100 deg^ cosmic shear dataset. 
From left: GaBoDS, RCS, CFHTLS, VIRMOS-DESCART, and 
combined. Each lA model is shown separately for each survey, 
plotted at the median redshift of each survey. From left: No lAs 
(horizontal crosses), linear alignment (circles), NLA model (diag- 
onal crosses) and halo model (triangles). lA points are offset from 
the median redshift for clarity. The combined constraints for the 
whole 100 deg^ dataset are shown in the panel on the right. This 
analysis varies erg leaving other parameters fixed at their fiducial 
values, including = 0.3. 
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Figure 4. 68% Confidence limits in erg - parameter space 
combining the four surveys in the 100 deg^ cosmic shear dataset 
for each lA model. Constraints where lAs are ignored in the mod- 
elling are shown by the dark, solid (blue) contour; constraints as- 
suming the linear alignment model are shown by the dark dotted 
(black) contour; the NLA model by light dotted (red) and the 
halo model by the light solid (green) contour. lA parameters are 
set to j4 = 1, /3 = 0, /r = 1. AH other cosmological parameters 
were fixed at their fiducial values. 



diction meaning that constraints using an lA model favour a 
larger value of erg to compensate for the reduced amplitude. 
The effect is more pronounced for the NLA model (diagonal 
crosses) than for the linear alignment model (circles) because 
suppression of the correlation function due to the GI term 
is stronger in the NLA. When the halo model is used (trian- 
gles) the best fit erg is larger than that when I As are ignored, 
but the difference can be smaller than for the NLA or LA 
models due to the contribution from the II correlations at in- 
termediate separations which partially cancel the GI terms. 
We note that the difference between the halo model con- 
straint and that when I As are ignored is least for GaBoDS 
and VIRMOS-DESCART. This is consistent with their rel- 
atively large number of low-redshift galaxies. This produces 
an increased II contribution at intermediate scales, lowering 
constrained as values. The effect is biggest for VIRMOS- 
DESCART because it has the greatest constraining power 
at these intermediate scales. 

Although lAs are an important cosmic shear systematic 
and a clear biasing is apparent, it is encouraging to note that 
the effect is not catastrophic: the best-fit values for almost 
all surveys and lA models lie within the 68% confidence lim- 
its of the results in which lAs are ignored. Because the un- 
certainties on as are smaller for the joint constraints (right 
hand panel of Fig. [3}, the total biases are now similar to 
the 68 per cent confidence limits. We use a more sophisti- 
cated model for how lAs affect cosmic shear constraints in 
Section [S] 

3.4 Impact on constraints on as and Qn-i 

We next calculate constraints from the cosmic shear data 
allowing both the fluctuation amplitude erg and the mat- 
ter density Q,n-i to vary. Other cosmological parameters are 



held at their fiducial values. The results are shown in Fig. 
|4] as 68% confidence contours in crg-nm parameter space. 
All results for the various lA models show the usual cos- 
mic shear degeneracy between erg and flm- As expected, at 
a fixed Slni — 0.3 we can read off the erg constraints shown 
in the previous figure. However the relative pattern between 
models depends on the assumed value of Qui and further- 
more different flm ranges are preferred by different models. 

The most interesting result is that from the halo model 
of I As (light solid (green) contour), which favours the high 
erg and low fim end of the degeneracy. Very roughly speak- 
ing the shear-shear correlation function data provide two 
pieces of information: an amplitude and a slope. The am- 
plitude constrains a degenerate combination of e7g and 57m 
but the slope of the correlation function partly breaks this 
degeneracy and determines the preferred range range of Sim 
along the degeneracy line. The most constraining data set is 
CFHTLS, that has tightest constraints on large scales. The 
effect of the halo model of lAs is to steepen the correlation 
function on these scales due to the contribution from the 
II one-halo term at intermediate scales (Note that this one- 
halo II term does not vary with cosmological parameters in 
our implementation so this conclusion stands at all positions 
in cosmological parameter space.) An increase in f2m flattens 
the matter power spectrum on the scales of interest, thereby 
increasing the amount of power on small scales relative to 
large scales. This means that an increase in Q^n leads to a 
more steeply decreasing correlation function. Thus to repro- 
duce the observed data in the presence of lA halo model 
contributions ilm must decrease, leading to the position of 
the halo model contours. 
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Figure 5. Comparison of the shear-position correlation function, Wg^-, from the lHirata et al] l|2007l ) data and the halo model prediction 
using fiducial cosmology, A = 1, /3 = 1.44, Lg = L4. Predicted correlation functions are shown for the one-halo term only (dotted), 
two-halo term only (dashed) and the total (solid). Blue galaxy data for each luminosity bin is plotted as crosses, red galaxies as diamonds. 
The shaded region shows distance scales excluded from our analysis. 



4 SHEAR-POSITION CORRELATIONS 



4.2 Correlation Functions 



In this Section we introduce a new data set which was pro- 
duced with the aim of measuring the lA signal. We use the 
galaxy shear - galaxy position two-point function that quan- 
tifies whether the shape of one galaxy points at the position 
of another galaxy. If we have a good understanding of how 
the galaxy positions are related to the mass distribution then 
this effectively constrains the lA power spectrum Pg . We 
compute predictions for the two point function using the 
halo model for lAs, then compare it to the data to measure 
lA model parameters. 



4.1 Data 

We use data measured from the SDSS galaxy survey, a 5- 
band (ugriz) photometric survey over tt sterad ians of the 
sky, w ith spectroscopic follow-up of 10^ objects jYork et al.l 
We use the galaxy position - g alaxy shear correla- 
tion fu nctions from the d ata described in Mandelbaum et al] 
(|2006l ) and calculated in lHirata et al.1 l|2007l ). These are cal- 
culated from the main spectroscopic sample (z — 0.05 - 0.2). 

This is split by luminosity and colour and each of these 
subsamples is used to trace the intrinsic shear field. The 
luminosity subsamples span L3 (one magnitude fainter than 
L*) to L6 (two magnitudes brighter). Each of these is divided 
into "red" and "blue" colour subsamples, defined using an 
empirically determined, redshift dependent separator, u — 
r — 2.1 -I- 4.22 in the observer frame. The full sample is 
used to tra ce the density field. E Uipticity measurements were 
taken from iMandelbaum et al.l (|2005.1 . 



The correlation function is calculated by taking pairs of 
galaxies where one galaxy is from the shear catalogue (lu- 
minosity and colour bin) and the other is from the den- 
sity catalogue (which is a composite of all the luminos- 
ity and colour bins). The ellipticity of the galaxy from the 
shear catalogue is calculated in a coordinate system aligned 
with the vector between the two galaxies in the pair. If 
the galaxy in the shear catalogue is aligned along the vec- 
tor joining the pair then this number will be positive. If 
it points perpendicular to this line then the ellipticity in 
this coordinate system will be negative. These ellipticities 
are averaged over all pairs, as a function of separation on 
the sky Vp and the separation along the line of sight 11, 
taking into account edge effec ts using a generalized Landy- 
Szalay l|Landv fc Szalavlll993h estimator. This is integrated 
over the range —60 < II//i~^Mpc < 60. We use the covari- 
ance matrices between bins in angular separation which were 
provided. 

We predict the shear-position correlation function Wg+ 



usmg 



ws+^—— I Pg;^i{k)J2(krp) kdk 



(17) 



(|Hirata et al.ll2007l . Eq. 23) which we assume is related to 
the galaxy correlation function by the galaxy bias bg 



^9+ ~ bgWS + 



(18) 



whe re we use ga laxy bias estimates from Table 3 of 
,Hirata et al.l (|2007l ) and we allow for the variation of galaxy 
bias with as as they describe. We use the generalised halo 
model as given in Eq. [15] and Eq. [16] and take values for 
the luminosity of each SDSS shear-position correlation func- 
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Figure 6. 68% confidence limits on A and /3, the amplitude and 
luminosity depe ndence of the lA po wer spectra. Constraints are 
plotted from the lHirata et all ||2007| ) shear-position data (dotted 
contour), the 100 deg^ shear-shear data (dashed contour) and 
joint constraints (solid contour), trg was marginalised over and all 
other cosmological parameters were fixed at their fiducial values. 



ti on bin from the min imum magnitude in each bin as given 
in iHirata et all |20o3) Table 1. We set Lo to be the lumi- 
nosity of the L4 bin. We set the red fraction fr to be unity 
for the red samples and zero for the blue samples, so that 
effectively we do not us e the b lue samples in the constraints. 
Following iHirata et al. 

I l|2007l ), we use only the data points 
for separations greater than 4.7h~^ Mpc where the bias can 
reasonably be assumed to be linear. The overall amplitude 
scaling parameter A and luminosity power law j3 are the two 
free parameters we consider. 

Fig. [5] shows the shear-position correlati on function, 
Wq+, f or each luminosity and colour bin in the lHirata et al.l 
l|2007l) data (crosses for blue galaxies, diamonds for red 
galaxies) . 

Also shown are predictions using the halo model, using 
for illustration A = fr = 1 and /? = 1.44, the best-fit /3 
valu e provided by the third row of the lower half of table 
4 m IHirata et alT(|2007^ ■ Separate correlation functions are 
plotted for the one-halo contribution (dotted line) , two-halo 
contribution (dashed line) and the total (solid line). As we 
would expect, the two-halo term provides the bulk of the GI 
signal at high separation, and falls ofT at small separation. 
An interesting feature of the one-halo correlation function is 
that as the separation is decreased it rises from zero, then 
falls to negative values, before rising s harply at low sepa- 
ration . This behaviour was discussed in lSchneider fc Bridle! 
l|2010l ) with reference to their Fig. 2 in their Section 2A. 
The grey-shaded regions on Fig. [S] show the scales which 
have been excluded from our ana l ysis o f the lA data, follow- 
ing the practice of IHirata et al. I (j2003). We see that essen- 
tially only the two-halo term is significant in the fits to this 
dataset. 



4.3 Constraints on Intrinsic Alignment 
Parameters 

In this section we use the shear-position correlation function 
data to constrain A and /?, the free parameters in the lA 



model. The dotted line in Fig. [6] shows the resulting 68% 
confidence limits. Both shear-shear and shear-position data 
are marginalised over erg. We chose to pivot the luminosity 
power law at the SDSS L4 luminosity bin (Lo corresponds 
to the magnitude of the L4 bin). The red galaxies are either 
in luminosity bin L4 or higher luminosity bins, therefore an 
increase in /3 will mostly increase the predicted lA signal 
relative to the data. The same is true for A, therefore we 
expect them to be anti-correlated, as seen in Fig. [S] We 
discuss the remaining lines on this figure in the following 
section, but note here that the shear-position correlation 
data is by far the most important dataset when constraining 
these two parameters alone. 

Our results agree well with those of lHirata et al.l (|2007l ) 
who find P = lAAtl ll when using aU the SDSS Main 
data (blue and red) simultaneously, with the same mini- 
mum separation in the correlation function. By construc- 
tion the best fit A is around unity, since the vast major- 
ity of the constraint comes from the two-halo term of the 
halo model which corresponds to the line ar theory predic- 
tion e stimated from Hirata fc SeliakI (|2004l ) which was fitted 
to the lHevman s ct al.' (2004*) data which gave a similar am- 
plitude as the [Mandolbaum ot al, (200Q) result. 

To allow a comparison between the cosmic shear data 
and the lA model used for the galaxy survey data values of 
fr and L were calculated as a function of redshift for each 
component survey of the 100 deg^ dataset using a condi- 
tional luminosity function model based on ICooravl (|2006l ) 
and outlined in appendix 1X1 

In the previous section we computed cosmological con- 
straints from cosmic shear data with and without consider- 
ation of lAs. The "No lA" results correspond to ^ = on 
this Figure, that is clearly ruled out by the SDSS data. The 
"halo model" results corresponded to A — 1, P — 0, fr — 1 
which is also ruled out by this figure but would be appro- 
priate if all the galaxies in the cosmic shear survey had the 
same luminosity and colour as the SDSS Red L4 sample. In 
the following Section we present a full joint analysis of both 
samples. 



5 JOINT SHEAR-SHEAR AND 

SHEAR-POSITION CONSTRAINTS 

In this Section we compare constraints on parameters from 
shear-shear and shear-position correlation functions and 
combine the constraints to produce limits on both lA and 
cosmological parameters. To make our previous analysis of 
the cosmic shear data more general it would be good prac- 
tice to allow the lA free parameters A and /3 to vary, then 
marginalise over them rather than assuming they were fixed 
to the basic halo model prescription. This would inevitably 
lead to a loss of constraining power. However, this loss can 
be made less severe by calculating joint constraints on lA 
parameters from cosmic shear and galaxy survey data. With 
this approach, the extra information about lAs provided by 
galaxy surveys is used to produce less biased cosmological 
constraints from cosmic shear. 
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Figure 7. 68% confidence limits on erg and A, the lA power 
spe ctrum amplitude p arameter. Constraints are plotted from 
the iHirata et al.l 1 120071) shear-posit i on co rrelation function data 
(dotted lines), the iBeni amin et al. 1 1I2OO7I') 100 deg-^ shear-shear 
correlation data (dashed contour) and joint constraints (solid con- 
tour). The lA luminosity dependence power law slope was fixed 
/3 = 1.44, and all other cosmological parameters were fixed at 
their fiducial values. 



5.1 Joint Constraints on Intrinsic Alignment 
Parameters 

We combine constraints from the two datasets by simply 
adding together their values. This is acceptable given 
that the surveys do not overlap significantly. We calculate 
values from both datasets as a function of both lA pa- 
rameters and cosmological parameters. For the shear-shear 
dataset the lA parameters enter via the lA power spectra 
(Eg. 1 151 and Ea. ll6p which are projected onto the sky (Eq. [S] 
and Eq. [6]) and added to the cosmic shear contribution to 
produce the full shear-shear power spectrum (Eq. [4]) and 
the correlation function is calculated by Eq. [S] The shear- 
position correlation functions are computed from the shear- 
density correlation function (Eg. llSp which is a inverse Han- 
kel transform (Eq. I17|) of the r elevant lA p o wer s pectrum 
(Eq.[6]). The bias values given in lHirata et al.l l|2007l ) are cal- 
culated from measurements of galaxy clustering (position- 
position correlation functions) and therefore depend on the 
assumed value of as- We take this into account to obtain a 
bias values for each erg value considered. As discussed at the 
end of Section r2.2l we fix as in the lA power spectra and al- 
low Q,rn to vary in the linear theory matter power spectrum 
contribution to the two-halo term. 

First we allow only erg and A, the amplitude of the lA 
signal, to vary, with the rest of the cosmological parameters 
set to their fiducial values and a fixed luminosity depen- 
dence power law slope /3 = 1.44 (the best-fit value given in 
iHirata et all [2007'). The results are shown in Fig. [7] as 68% 
confidence contours in as-A parameter space. Contours are 
shown for the shear-position correlation functions (nearly 
vertical lines) shear-shear correlation functions (nearly hor- 
izontally elongated contour) and the combined constraint 
(roughly at the intersection). The shear-position correlation 
function constraints on the lA amplitude parameter A for 
the fiducial ag value are as expected from Fig. [6] for the fixed 
P value, with A significantly larger than zero. In the analy- 



sis of the shear-position correlation function data ag is held 
fixed in in the calculation of all our lA models. However, 
the variation of galaxy bias as a function of as produces 
the expected degeneracy with A. Increasing as above the 
fiducial value decreases galaxy bias, requiring a greater A to 
compensate and vice-versa. 

The shear-shear correlation functions do themselves 
place some constraint on the lA parameter A, preferring a 
range of order unity. A negative value of A would correspond 
to galaxies pointing in the opposite direction to that used 
in the standard models. The linear alignment model (two- 
halo term) would contain galaxies pointing perpendicular 
to the tidal stretching expected by the gravitational poten- 
tial curvature. The one-halo picture would contain galaxies 
which are aligned tangentially to the center of the halo. We 
described earlier a rough picture of shear-shear correlation 
function constraints in which essentially the data measure 
the amplitude and slope of the correlation function. The 
amplitude essentially fixes a degenerate combination of erg 
and A, and the constraint on A must therefore come from 
the shape of the correlation function. For the fiducial ilm, 
Fig. [7]tells us that a large lA contribution to the shear-shear 
power spectra distorts them too much. This can be under- 
stood by examining Fig. [21 in which the data points fit well 
with the shape of the lensing-only ("No lA") predictions, 
whereas the halo model predictions tend to be more curved 
over the scales probed. 

The direction of degeneracy between A and as from 
shear-shear information alone can also be understood in 
terms of the shear-shear correlation function datapoints in 
Fig. [2] In general there is a balance of effects between the 
II and GI contributions. If we imagine a universe in which 
the GI term did not exist but the II term did, an increase in 
the lA amplitude parameter A would add power to the pre- 
dicted shear-shear correlation function. To keep the predic- 
tions consistent with the data points it would be necessary 
to decrease erg to reduce the lensing contribution. This would 
give a negative correlation between A and ag , seen for larger 
positive A values. In this unphysical universe containing only 
II terms, the direction of degeneracy would appear reversed 
for negative A because A appears as a squared quantity in 
the II terms. Indeed we do see the direction reversed in this 
figure. In the physical Universe the effect of the GI contri- 
bution is to break the symmetry slightly, and prefer slightly 
more positive A values since it causes a slight cancellation 
in the lA effect for positive A values which fits better to the 
shape of the correlation functions that an addition of power. 

The great complementarity of the two datasets is most 
clearly illustrated in this figure. The main constraints on 
cosmology come from the shear-shear correlation data and 
the main constraints on lAs from the shear-position data. 
The joint constraints focus at the intersection between the 
two relatively degenerate constraints, as expected. 

The analysis was repeated, allowing A, (3 and ag to vary. 
We marginalised over ag to produce the results shown in 
Fig. [6] All other parameters were fixed at their fiducial val- 
ues. The shear-shear correlation function data (dashed) pro- 
vides much weaker constraints in the A-/3 plane- partly due 
to marginalisaion over ag in the power spectrum but primar- 
ily because the cosmic shear data is not broken into lumi- 
nosity bins, so has very little power to constrain /3. Even so 
we can distinguish the same degeneracy direction as present 
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Figure 8. 68% confidence limits on erg for each survey within the 
100 deg^ shear-shear correlation function dataset (from left: Ga- 
BoDS, RCS, CFHTLS, VIRMOS-DESCART and combined). For 
each survey constraints shown, from left: without lAs (crosses); 
using the halo model with A = I, P = 1.44 (circles); halo model, 
marginalised over A and /3 (stars); joint shear-shear and shear- 
position constraints, marginalised over A and /3 (diamonds). The 
right-side panel shows joint constraints from the four surveys. All 
other cosmological parameters were fixed at their fiducial values. 

in the shear-position data and positive A is favoured at 68% 
confidence. Most of the constraining power in the joint con- 
straints (solid) comes from the shear-position data (dotted). 

5.2 Joint Constraints on Cosmology 

Here we revisit the constraints on cosmological parameters 
from Section 13.21 showing how they change when we allow 
the lA parameters (A and /3) to vary and be marginalised 
over. Joint shear-shear and shear-position constraints are 
also presented to show constraints on cosmology when lAs 
are taken into account as well as possible given current mod- 
els and data. 

Fig. [8] shows constraints on the amplitude of matter 
fluctuations erg for each shear-shear survey separately as a 
function of survey median redshift, and also the joint con- 
straint for all shear-shear surveys in the right hand panel. 
The blue crosses in Fig. [8] are the same as those in Fig. [3] and 
show constraints on ag when lAs are ignored. In Fig. O we 
illustrated the effect of lAs on the different lA models using 
a simple version of the I A power spectra in which A = 1, 
(3 = and we assumed /r = 1. In Fig. [5] we use only the 
halo model and use the values for L and fr as a function of 
redshift as discussed in Appendix]^ For the circles in Fig. [8] 
we fixed A = 1, = 1.44, and we see that the as results us- 
ing the halo model are all now lower than before. The result 
of using the more realistic L and fr values instead of /3 = 0, 
/r = 1 is to lower all the lA contributions. This is to be 
expected since fr is always less than or equal to unity in 
the more realistic model, and in addition fr is particularly 
small at low redshifts where faint satellite galaxies dominate 
and the I A power spectra are largest. The biggest modifi- 
cation to the predicted correlation functions is the reduced 
GI contribution, especially on smaller scales. Therefore the 
total shear-shear correlation is increased and the fitted erg 
is reduced to compensate. The trends with different shear- 



shear surveys reflect those already discussed in reference to 
Fig. [3l As a result, the joint constraint (right panel of Fig. [S]) 
shows the halo model with A = 1, P = 1.44 (circle) is now 
slightly lower than when lAs are ignored. This is consistent 
with the contours shown in Fig. [7] which have a lower erg at 
A = l than A = 0. 

We next allow A and /3 to vary within the halo model, 
and we marginalise over both parameters with fiat wide pri- 
ors using only the shear-shear correlation data. The results 
are shown by the stars in Fig. [8] As expected, the uncertain- 
ties are larger than when these parameters are kept fixed. 
The best fit values are about the same as when no lAs are in- 
cluded. We would expect this from Fig. [7] in which the A = 
line cuts a cross-section through the banana-shaped degen- 
eracy contour from shear-shear correlations alone which is 
roughly in the middle of the ug range allowed by the whole 
contour. 

Joint constraints from the combined cosmic shear and 
lA data (triangles) gives tighter constraints on ag than from 
cosmic shear data alone, as expected. The constraint from all 
shear-shear surveys (right hand panel) is consistent with the 
contours in Fig. [71 where it is clear that the joint constraint 
(solid line) gives a very similar erg value to the cosmic shear 
constraint without lAs (a slice through the cosmic shear 
contour at A=0). 

The variation between different lA treatments are 
most strongly pronounced for the case of the VIRMOS- 
DESCART survey (highest median redshift of the four sur- 
vey points in Fig. [8]). Although it has the highest median 
redshift, we have seen from Fig. [l] that it also contains a 
large fraction of low-redshift galaxies and hence most phys- 
ically close galaxies. This increases the II lA signal, that 
is significant for the range of scales probed by VIRMOS- 
DESCART. This boosts the amplitude of the predicted 
shear-shear correlation function points, requiring lower erg 
to compensate. This behaviour produces the larger bias ob- 
served in VIRMOS-DESCART on going from no lA to the 
halo model. 

The new as value can be higher or lower than when lAs 
are ignored (horizontal crosses), depending on the dataset 
in question. However all the new cosmology constraints (tri- 
angles) lie easily within the 68% confidence limits of the 
constraint without lA, including the constraint for all shear- 
shear surveys combined. 

We repeated the analysis for constraints in erg-Sim space. 
All other cosmological parameters were held at their fiducial 
values. Fig. [9] shows the 68% confidence contours for cosmic 
shear data without lAs (blue solid contour), that are the 
same as the dark solid lines in Fig. [4] The halo model with 
fixed lA parameters (green dotted contour) shows similar 
behaviour to that discussed in Fig. IH favouring the high 
erg and low Qm region. Marginalising over the lA ampli- 
tude and luminosity-dependence parameters, A and /3 (red 
dotted contour), relaxes the constraint, extending the de- 
generacy further into the high erg, low Qm region. The joint 
constraints (black solid contour) using both shear-shear and 
shear-position information show a preference for low as and 
high Qm but are still consistent with our fiducial cosmol- 
ogy. Th e preference for high er Jim values occurs because, as 
noted in lHirata et al.l ()2007l ). the shear-position data favours 
slightly steeper correlation functions than predicted by a 
standard fiducial cosmology (see also the L6 panel of Fig. [5]). 
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Figure 9. Combined 68% confidence limits in crg-Qni parameter 
space for tfie 100 deg^ siiear-sfiear correlation function dataset. 
Constraints stiown without I As (blue solid contour); using the 
halo model with A = 1, 13 = 1.44 (green dotted contour); halo 
model, marginalised over A and /3 (red dotted contour); joint 
shear-shear and shear-position constraints, marginalised over A 
and /3 (black solid contour). All other cosmological parameters 
were fixed at their fiducial values. 

Increasing Qm produces a correlation function that falls off 
faster at larger scales because the linear theory matter power 
spectrum peak is moved to the right (smaller scales) on in- 
creasing Qn-i due to the earlier epoch of matter-radiation 
equality. 

This new constraint on cosmology from cosmic shear 
self-consistently takes into account lA effects and we find 
fim = 0.18 ± 0.09, (78 = 0.88 ± 0.18 at 68 % confidence, as 
opposed to the results in which lAs are ignored, that give 
fim = 0.20 ± 0.10, o-g = 0.86 ± 0.21. When we constrain 
simultaneously with shear-shear and shear position data we 
find r^m = 0.29 ± 0.09, erg = 0.74 ± 0.13, after marginalising 
over lA parameters. There is a strong degeneracy between erg 
and Qn-i when using current cosmic shear data alone, so it is 
better to look at the full two-dimensional contours (Fig. ^ 
or at cross-sections. When we fix ilm = 0.3 our best self- 
consistent shear-shear result including lAs (right hand panel 
of Fig. [8| is ug = 0.75 ± 0.04 as compared to that when lAs 
are ignored of ug — 0.75 ± 0.03. A joint constraint from 
shear-shear and shear-position data, marginalising over lA 
parameters, gives erg = 0.75 ± 0.03. 



6 CONCLUSIONS 

lAs are expected to be an important source of systematic er- 
ror in cosmic shear measurements if completely ignored. Sev- 
eral previous constraints on cosmology from shear-shear cor- 
relation functions have considered the possible levels of con- 
tamination from lAs, but the recent constraints presented 
have not included the impact of lAs (iBeniamin et al.ll2007l : 
et al. I I2OO8I : ISchrabback et aLlbOoi T These surveys are 
therefore referred to as "cosmic shear" surveys. 

In this paper we use three physically motivated lA mod- 
els, of increasing complexity, to demonstrate how constraints 
on the amplitude of fluctuations erg and matter density Qm 
axe changed depending on the assumptions about lAs. The 



size of the change in erg depends on the redshift distribu- 
tion of galaxies in the shear-shear survey and on the range 
of scales probed. In particular this affects the interplay be- 
tween the two different types of lA (II and GI) which have 
competing effects on cosmological constraints. A simple ex- 
amination of lAs which assumes all galaxies are like the 
SDSS Main Red L4 sample suggests that erg has previously 
been underestimated by about one standard deviation. The 
underestimate appears to be larger using the non-linear mat- 
ter power spectrum in the linear alignment model for lAs 
than when using the linear matter power spectrum. On us- 
ing the halo model for lAs the bias depends significantly on 
the range of scales probed by the survey due to the larger 
effect of II at intermediate scales which biases the measured 
(jg downwards when taken into account. 

We have performed the first constraints on multiple lA 
model parameters from observed shear-position correlation 
functions using a physically motivated model for lAs. We 
also show the first simultaneous constraints on lA and cos- 
mological parameters from either shear-position or shear- 
shear correlation functions. As expected, the majority of the 
constraint on cosmology comes from the shear-shear corre- 
lation functions, and the majority of the constraint on lAs 
comes from the shear-position correlation functions. 

We have used a motivated model for luminosity and 
colour evolution of shear-shear galaxy samples as a func- 
tion of redshift to improve our constraints on cosmology 
from shear-shear correlation data in the presence of the halo 
model for I As. In general this reduces the overall effect of I As 
on cosmology, which we attribute to the decreased contribu- 
tion from lA at low redshifts. We consider constraints on the 
amplitude of fluctutions erg from shear-shear data alone with 
fixed lA amplitude and luminosity dependence parameters, 
and compare this with constraints after marginalising over 
these parameters. As expected the constraints are weaker 
after marginalisation, but are not biased signiflcantly. This 
is due to the relatively fiat degeneracy between erg and the 
lA amplitude arising from the competing effects of GI and 
II terms, and the ability of shear-shear information alone to 
place some constraint on the lA amplitude if other cosmo- 
logical parameters are held fixed. 

The model we used for the luminosity and colour evo- 
lution of the shear-shear galaxy sample could be improved 
in its complexity, for example by allowing a population of 
faint blue galaxies at high redshift. Furthermore we assumed 
the shear-shear galaxy sample had a single luminosity and 
colour at each redshift, whereas a more sophisticated anal- 
ysis would integrate over ranges in each. 

A joint analysis of shear-shear and shear-position cor- 
relation functions has been discussed as a promisi ng tool 
for removing lAs from future cosmic shear da tasets jZhanel 
l2008l : lBernsteinll2009l : |joachimi fc Bridlell2009 l'). We perform 
the first joint analysis of these correlation functions on ob- 
servational data. The uncertainties on the amplitude of fluc- 
tuations erg are now reduced to the size found when lAs are 
ignored. The estimated value of erg is either decreased or 
increased relative to the case where lAs are ignored. How- 
ever, f or the combined constraints from the Benjam in et al.l 
(|2007l ) shear-shear data compilation the overall erg value is 
barely changed on performing the full joint analysis of shear- 
shear and shear-position correlation functions with lAs, rel- 
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ative to the conventional analysis in which lAs are ignored 
and only shear-shear information is used. 

The constraints in the two-dimensional parameter space 
of erg and Sim show a bigger change, that we have attributed 
to the slope of the shear-position correlation function pre- 
ferring higher f2m values. This raises the question of how to 
incorporate cosmological parameters into lA predictions. For 
example we could have considered the lA power spectra from 
our fiducial cosmological model as template functions which 
are allowed to vary only in amplitude. It also places a new 
burden of accuracy on measurements of the shear-position 
correlation function, that in turn depend on measurements 
of galaxy biasing for their interpretation. 

We have not repeated the full analysis of 
iBeni amin et"al] (|2007t ). who consider several other potential 
sources of systematic uncertainty including photometric 
redshift and shear calibration errors. We defer this to 
future investigations that could incorporate the very latest 
cosmic shear data. The recent re-anal ysis of the COSMOS 
cosmic shear data (|Massev et al ] |2007l ) bv lSchrabback et al.l 
l|2009D investigates a number of potential cosmic shear 
systematics including I As. The photometric redshifts are 
used to produce the first cosmic shear tomographic cross- 
correlation functions which allow a much better control of 
lAs. Furthermore they investigate the effect of removing red 
galaxies from the cosmic shear sample which should help 
to reduce the effect of lAs further. This would have to be 
carefully modelled in the predictions of the lA contribution 
to the shear-shear correlation functions. 

We believe that future cosmic shear analyses should in- 
clude a model for lAs and use the existing galaxy survey 
data to help constrain this simultaneously with cosmology, 
as set out in this paper. The modelling in this paper could be 
improved by extending the flexibility in the lA model, for 
example by allowing the one- and two-halo terms to vary 
independently in amplitude, and allowing an unknown am- 
plitude for the lAs for blue galaxies. Furthermore the ampli- 
tude parameters could be free functions of redshift and even 
scale. The lA models we have used specify simultaneously 
the two types of lA, II and GI. However, a more general 
model of lAs can lead to independe nt uncertainties in each 
of the source terms of these effects (|Bernsteinll2009l . see the 
bias and correlation terms of the intrinsic convergence with 
the lensing potential in). 

However, if such a large degree of flexibility were al- 
lowed with current shear-position information then the con- 
straining power of current cosmic shear data would be sig- 
nificantly reduced. Therefore it is important to obtain bet- 
ter measurements of the shear-position correlation functions, 
and/or improve models of lAs to reduce the number of free 
parameters. In this analysis we have used measurements of 
the galaxy shear-position correlation function made from a 
spectroscopic galaxy redshift survey. It should be possible 
to make these measurements from the same data used for 
the cosmic shear survey itself, however this has not yet been 
done. This will require a careful understanding of the survey 
mask as well as a measurement of the galaxy biasing to al- 
low the shear-position correlation function to be related to 
the shear-density correlation function. In addition there are 
constraints in the literature from the shear-shear correlation 
functions of low redshift galaxy samples as well as for higher 



redshift spectroscopic galaxy samples which can be used to 
provide further constraints on the lA contributions. 
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APPENDIX A: MODELS FOR MEAN 
LUMINOSITY AND RED GALAXY FRACTION 

To match the luminosity and colour dependence of the lA 
model in the cosmic shear data we use a model for the con- 
dition al luminosity function (CLF) closely following ICooravl 
(|2006l ). The conditional luminosity function $(L[M, z) gives 
the mean number of galaxies in a dark matter halo of mass 
Al at redshift z with luminosity in the range [L,L + dL]. 
This function provides a way to add luminosity dependence 
to the halo model description of galaxy statistics. 

Following the usual practice in the halo model, the CLF 
model is constructed from separate models for the central 
and satellite galaxies. The first step is to specify the mean lu- 
minosity of the central galaxy in a halo of mass Al, Lc{AI, z), 
which we model as a brok en power-law in mass according 
to Eq. 3 of ICooravl (|2006l ): where the mean luminosity of 
the central galaxy increases quickly with mass for low halo 
masses and then fiattens off for larger halo masses. The total 
luminosity of a halo should include contributions from satel- 
lite galaxies for masses greater than some threshold, Alsat, 



Ltot(M, z) = 



Lc{Ai,z) M=^M,at(z) 

(Al) 

where Ps = 0.8, A'h^t{z) = Msaterf(5z) and Msat = 
10^^ h~^AdQ. This particular z dependence is an ad-hoc ad- 
justment to allow lower mass halos to host satellites at low 
z. 

In the simplest model the CLF for cen tral galaxies is 
just a delta-function centered at Lc, although lCooravl (|2006l ) 
suggests a "rounded" step to allow for scatter in the mini- 
mum mass halo hosting a galaxy, 

^"""(LIM, z) = Bh (M - Mcen-cut) Log-Normal(L | , (Tc), 

(A2) 

where Oc = 1.15, Qh is the Heaviside step function and 

Afccn-cut = lO^/i-^MQ. 

The satellite CLF is modeled as a power-law in lumi- 
nosity with maximum luminosity equal to half the central 
galaxy luminosity (so that satellites do not dominate the 
total luminosity of the halo), 

c&=^*(L|Af, z) ^Qh(^- I^j A{M, z) J. (A3) 

The amplitude, A{A4,z), is determined by requiring 
J^'""'"' dLL^'"'\L\M, z) = Ltot - Lc, giving. 



A{M,z) 



Ltot{A'I,z) - Lc{M,z) 
I {Luiax, Lc{M, z)/2) — L 

min ) 



(A4) 
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Figure Al. Mean luminosity model versus redshift for each of 
the 4 lensing surveys. 



For modeling the luminosity dependence in the 100 sq. deg. 
lensing data sets, we set the integration limits Lmin,imax 
to match the app a rent m agnitude limits given in table 1 
of iBeniamin et aP (l2007l ). Note the luminosity integration 
limits become redshift dependent when matching fixed ap- 
parent magnitude limits. For the upper integration limit we 
further set Lniax(M, z) = min {Lc{M, z)/2, Lmax,survoy (z))- 

The HOD functions giving the number of central and 
satellite galaxies in a halo of mass M are obtained by inte- 
grating the CLFs over luminosity, 



dL 



\L\M, 



(A5) 



Integrating over halo mass, we can relate the CLFs to 
the redshift distribution of the galaxy sample as. 



dM n{M, z) 
X (iVrXAf,z) + iVr(M,^)), 



(A6) 



where is the comoving distance as a function 

of redshiflu. We ma tch the redshift distributions from 
IBeniamin et al.l (|2007l ') with the CLF prediction to fit the 
(redshift-dependent) normalization of the central galaxy lu- 
minosity Lc(M, z). 

With the parameters of our CLF model matched to the 
magnitude limits and redshift distribution of each lensing 
survey, we calculate the mean luminosity of each galaxy sam- 
ple as. 



{L) {z) 



J dM n(M, z) J dLL^{L\M,z 
J dM n{M, z) J dL $(L|M, z) 



(A7) 



We show our model for (L) (z) for each of the four lensing 

surveys in Fig. lAll 

We again follow [Cooravl (|2006l l to model the fraction of 



^ assuming zero spatial curvature. 



Figure A2. Red galaxy fraction versus redshift for each of the 4 
lensing surveys. The top panels show the red fraction for central 
and satellites separately in our halo model while the bottom panel 
shows the red fraction for all galaxies. 



red central galaxies as a function of halo mass as a "rounded" 
step function, 



1 



/red, con (M, z) = ficd,ccn{z)- |^ H- Crf 



^ ln(M) -ln(Mecn) 

^rcd,ccn 



(A8) 

This approaches the constant value /red, con (z) as M becomes 
large and goes to zero for small M. Note this has no lumi- 
nosity dependence. 

The fraction of red satellites in a halo of mass M is 
modeled as, 

/rcd,sat(Af, L, z) = 3rod,sat(^) {g{M,z) /l(L, z) ) + /rcd,sat (z) 

(A9) 

where g{M, z) and h{L, z) are "rounded step functions" in 
mass and luminosity, respectively. For low masses and lumi- 
nosities, the red satellite fraction is the constant /red, sat (z) 
while for high masses and luminosities it is the different con- 
stant 2grcd,sat(z) -f /rcd,sat(z). For simplicity, we set /red,con, 

/red, sat, and (?red,sat to the z-indepeudeut constants 0.6, 0.4, 
and 0.2. 

From these pieces, the fraction of all red galaxies as a 
function of redshift is. 



/rod(z:) = 



dMn{M,z) / dL 



N,{z) 

X [r'="(L|M,2)/red,con(M,2) 

+ 4''"'(L|M,2)/rod,sat(M,L,2)], 



(AlO) 



where Ng{z) is the comoving number density of galaxies in 
the sample. 

Our models for /red (2) for each of the 4 lensing surveys 
are shown in fig. IA2I Ignoring redshift evolution in /rod, con 
and /red, sat is almost certainly incorrect and, in combination 
with the fixed magnitude cutoffs in our model, tends to select 
excess red galaxies at high redshift (when normalized to give 
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reasonable red fractions at low redshift). Because the I A 
halo model includes contributions only from red galaxies, we 
consider this model for /red to give an upper bound on the 
possible lA contamination to the weak lensing data. This is 
suitable for the goal in this paper of assessing the impact of 
lAs as systematic errors in lensing surveys. To improve the 
predictive power of our lA halo model in the future it will 
be important to better understand the redshift evolution of 
the red fraction of galaxies. 



